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ABSTRACT 

Using the IRAM Plateau de Bure Interferometer, we have searched for CO (3-2) emission from the 
dusty Lyman break galaxy Westphal-MDll at z = 2.98. Our sensitive upper limit is surprisingly low 
relative to the system's 850 ^m flux density and implies a far-IR/CO luminosity ratio as elevated as 
those seen in local ultraluminous mergers. We conclude that the observed dust emission must originate 
in a compact structure radiating near its blackbody limit and that a relatively modest molecular gas 
reservoir must be fuelling an intense nuclear starburst (and/or deeply buried active nucleus) that may 
have been triggered by a major merger. In this regard, Westphal-MDll contrasts strikingly with the 
lensed Lyman break galaxy MS 1512-cB58, which is being observed apparently midway through an 
extended episode of more quiescent disk star formation. 

Subject headings: cosmology: observations — galaxies: ISM — galaxies: starburst 



1. INTRODUCTION 

Among z > 2 source populations, Lyman break 
galaxies (LBGs) dominate in sheer numbers of objects 

■ with spect roscopic redshifts ( Stcid el et al.l 119991 1200,'^ 
, lLehnert~ Bremer 2003'). As a result, they define a nat- 
ural test set for galaxy evolution models (e.g., Baugh 

; et al. 1998; Mo, Mao, & White 1999; Somerville, Pri- 
mack, & Faber 2001). Different authors have reached 
different conclusions about the nature of star forma- 
tion in LBGs. The population's stron g bias relative to 
the u nderlying dark mat ter fluc tuations llGiava lisco et al. 

■ 119981 lAdelbereer et alJ 119981: lArnouts et all 19991. at 
a level apparently correlated with rest-UV luminosity 

: HGiavaliscofc DickinsonlEioTt iFouca.ud et Hni2()()^. a,r- 
gues that LBGs are forming spheroids at the centers 

' of the most massive z ~ 3 halos. In contrast, a 
high close pair fraction suggests that the most mas- 
sive halos host more than one LBG (Wechsler et al. 

' 2001; Bullock, Wechsler, & Somerville 2002), support- 
ing scenarios in which LBG starbursts are brief and 

' merger-driven { Lowcnthal et al. 1997; Ko latt et alJll999t 
iSomerville et abi200Ll . Because the star formation his- 
tories in these two pictures differ starkly, it is of great in- 
terest to establish observationally how LBGs' star forma- 
tion rates, stellar population ages, stellar and gas masses, 
and dynamical timescales relate to each other. In partic- 
ular, fully understanding star formation in LBGs requires 
characterizing the molecular gas from which their stars 
form; this motivates searches for CO line emission in at 
least a few representative systems. 

The logical first target for such an effort was the gravi- 
tationally lensed system MS 1512-cB58 (hereafter cB58), 
the brightest LBG known. Baker et al. (2004; B04) suc- 
ceeded in detecting cB58's CO(3-2) emission after prob- 
ing to greater depths than were reached by two previously 
published nondetections. As a second target, we chose 

^ Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by the INSU/CNRS 
(France), the MPG (Germany), and the IGN (Spain). 

^ Also at Department of Physics, 366 Le Conte Hall, University 
of California, Berkeley, CA 94720-7300. 



Westphal-MDll^ (hereafter WMDll), the dustiest LBG 
known. This z = 2.98 system was the strongest 850 /im 
source in the Submillimet er Common-User Bolom eter 
Array (SCUBA) survey of iChaoman et all pOO0|) and 
was revealed in high-resolution Hubble Space Telescope 
(HST) and Keck imaging to be a composite system in- 
cluding one Extremely Red Object (ERO) component, 
two bluer ones, and a web of diffuse blue emission 
l|Chanman et al.ll2002D . Sin ce WMDll 's 5's.';n resem- 
bled that observed for cB58 llvan der Werf et all llOOl). 
we reasoned that the two galaxies' CO emission lines 
should also be comparably strong. Instead, we report 
an upper limit on WMDll's CO(3-2) emission based on 
sensitive millimeter interferometry. This paper assumes 
Hq = 70/10.7 km s-iMpc-\ f^Af = 0.3, and 17a = 0.7 
throughout. 

2. DATA ACQUISITION AND REDUCTION 

We observed WMDll with the IRAM Plateau de Bure 
Interferometer (PdBI: Guilloteau et al. 1992) from 2003 
April through September. The array included four to six 
15 m diameter antennas; during our observations, these 
were arranged in variations on a compact D configura- 
tion with baselines ranging from 24 to 113 m. The 3 mm 
receiver on each antenna was tuned in single-sideband 
mode, giving typical system temperatures of 100-200 K 
in the lower sideband relative to the reference. To ob- 
serve WMDll, we adopted as a pointing center the posi- 
tion (aj2(Mi^I4H8™09J31and (5j2ooo = -h52°22'01.3") 
hsted bv isteidel etTTI (l2f)f)l a^nd use d for the SCUBA 
observations of 'Chaoman et al.' (2000). For the zhii = 
2.9816 reported by Pettini et al. (2001), the CO(3-2) ro- 
tational transition is redshifted to 86.85 GHz. We de- 
ployed four correlator modules at this frequency, giving 
a total of 560 MHz (1933 km s~^) of contiguous band- 
width at 2.5 MHz (8.6kms~^) resolution. Although we 
also tuned each antenna's 1 mm receiver for simultaneous 
242 GHz observations of thermal dust emission, the sum- 
mer weather conditions caused prohibitively high phase 

^ This designation supersedes the earl ier Westphal-MMDll in 
the master catalog of lSteidel et alj 120031) . 
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noise and prevented us from using these data. 

We calibrated t he data using the CLIC ro utines in the 
GILD AS package ijGuilloteau fc Lucasl200Cl|) . Phase and 
ampUtude variations within each track were removed by 
interleaving observations of the quasar J1419+543 (2.0° 
from WMDll on the sky) every 30 minutes. Passband 
calibration used one of several bright quasars. The flux 
scale for each epoch was set by comparing J1419+543 
with the passband calibrators and model sources CRL618 
and MWC349, whose flux densities are regularly moni- 
tored at the PdBI and IRAM 30 m. From the variance 
of multiple measurements, we estimate the accuracy of 
our final flux scale to be ~ 15%. Before constructing 
uv tables from our calibrated data, we smoothed them 
to 5 MHz (i.e., 17.2 kms~^) resolution. After editing for 
quality, we were left with the (on-source, six-telescope 
array) equivalent of 42 hr of data. We inverted the data 
cube without deconvolution using the IMAGR task in the 
NRAO AIPS package (van Moorscl. Kcmball. & GreisenI 
1199(f). With natural weighting, the synthesized beam is 
8.6" X 5.9" at P. A. 61°; the noise level per channel is 
0.7 mJy beam" ^ across most of the bandpass before ris- 
ing slightly at the high-velocity (low-frequency) end. 

After we had already acquired the PdBI data, we 
learned (N. Reddy, private communication) that a com- 
prehensive registration of radio, optical, and X-ray 
sources in the Westphal field implies a true position for 
WMDll some 2.5" northeast of its published coordinates 
(at ajaooo = 14''18"'09.78'* and Snooo = -h52°22'03.6"). 
This offset implies- pointing errors aside- that for the 
14.7" HPBW of S CUBA at 850 uin the Ssso = 5.5 ± 
1.4 mJy reported bv lChapman et all (2000) should trans- 
late to Ss5o = 5.9 ± 1.4 mJy at the actual position of the 
LBG. The offset also means that any CO (3-2) emission 
from WMDll should not lie at the PdBI phase center. 

With cursory inspection showing no strong source, we 
searched for line emission by convolving the 17.2 km s~^ 
cube with seven- and sixteen-channel boxcars. The 
resulting cubes have 121kms~^ (276 kms"^) chan- 
nels, each of which overlaps its nearest neighbors by 
104kms"^ (259 kms"^). The 121 k ms"^ width was cho- 
sen to match cthq = 53±5kms-i (jPettini et alJl200l|) 
by analogy with cB58, in which cthq — crco (B04). The 
276kms~^ width allows for a ratio ctco/chq — 2.2, the 
highest value seen for crco/fpaa among twelve Ultra- 
Luminous InfraRed G alaxi es (ULIRGs) ob s erved by both 
iSolomon et all l(l99l and iMurnhv et all l|200l .^ Nei- 
ther convolved cube shows a significant source at the po- 
sition of WMDll. Corrected for primary beam response, 
the 3cr upper limit of Fco(3-2) < 0.11 (0.17) Jykms"^ 
implies a line luminosity L'^o{3-2) — ^-^ (^-4) ^ 
10^ /i^^Kkms-ipc^ for Avqo = 121 (276) km s^^. 

3. WMDll CONTAINS A COMPACT DUST SOURCE 

In view of the encouraging intelligence available before 
we began our campaign, the failure to detect WMDll 
is rather surprising. Empirically, its line/continuum ra- 
tio (calculated here as -Fco(3-2)/'5'85o) of < 29kms~^ is 
< 1/3 of that seen for cB58 and < 1/7 of the mean 

This gives a conservative upper limit to the velocity width of 
CO(3-2) emission from WMDll: the data of Murohv et al. (2001) 
imply a mean (cco/crpaa) = 1.1, while the data of .Armus et all 
J1989I) for a smaller ULIRG sample give {o'co / '^Ha) = 0-9. 



for four z ~ 2.5 submillimet er galaxies (S MGs) with 
published CO(3- 2) detections l|Fraver et al.l lT998. 199i 
iNeri et al.lf2003j) . WMDll probably does not have a 
strong CO(3-2) line lying outside the bandpass of our 
PdBI observations. Its zhii is known to ±20kms~^ 
(M. Pettini, private communication) an d is nearly iden- 
tical for its blue and ERO components l)Chapman et alJ 
l2002j) . Relative to this zhii, its interstellar absorption 
lines are blueshifted by only 170 kms ~^ and its Lya emis - 
sion redshifted by only 230 km s'^ (jPettini et al.l l2001'l. 
Given such a limited velocity range (and the modest 
zhii — zco ^ 200kms~^ seen in cB58 by B04), it seems 
very unlikely that a molecular line could lie at such a 
high (> 966 km s"^) offset from WMDll's systemic ve- 
locity as to escape detection. 

An alternate explanation for our CO(3-2) nonde- 
tection would p osit t hat the dust emission seen by 
iChapman et"all l|2000() comes from a highly obscured 
companion or background gala xy rather than from 
WMDll itself. iSmail et all l|2002D show the surface den- 
sity of SMGs with 5'85o > 5 mJy to be 880^330 deg"^. 
might then expect the original survey of 16 LBG fields by 
IChapman et al.l l|2000(l to have turned up 0.3to;i sources 
in this brightness regime at random. While such odds 
are not negligib le, WMDll's unusually red TZ — color 
(|Shaplev et al.l 120011. together with the systematic de- 
pendence of 1.2 mm flux density on this color in a larger 
LBG sample (A. Baker et al. 2004, in preparation), sug- 
gests that some association of the 850 /xm source with the 
rest UV/optical counterpart is reasonable. 

Ruling out a nearby companion to WMDll takes more 
effort. If we no longer require that the 850 /Ltm emis- 
sion come from WMDll per se, we can look for CO 
sources that differ sharply from the LBG in position 
(e.g., Neri et al. 2003), velocity, and/or line width. We 
have therefore smoothed the data cube to several velocity 
widths (up to ^ 700 km s""'^) and searched for compan- 
ions within ~ 10" of the PdBI phase center (beyond this 
radius, a source would have to be unrealistically bright 
to have produced 5.5 mJy at the SCUBA pointing cen- 
ter). In each smoothed cube, the strongest candidate 
companion has < 3.8(J significance, lies near a bowl of 
comparable (negative) significance, and lacks an optical 
counterpart when the pipeline-reduced TJS'T imaging of 
IChapman et alJ l|2002|) is re-examined. In most cases, 
the implied Jco(3-2)/'S'850 is even lower than 29kms~^ 
owing to the larger correction that must be made to S'sso 
for a "source" farther from the phase center. We do not 
view any of these candidates as real. It remains possi- 
ble that a dusty companion lies close enough on the sky 
to WMDll to yield a strong 850 /im detection but far 
enough away in velocity to escape our PdBI bandpass. 
However, this scenario is dismayingly ad hoc, in that it 
divorces a high dust luminosity- at a level attained only 
by galaxy mergers in the local universe- from a system 
with a clear merger morphology. 

It remains to explain how WMDll can generate a large 
dust luminosity from a small reservoir of molecular gas. 
Exactly how much of a challenge this is depends on the 
value of the galaxy's IR luminosity, which we can calcu- 
late from its 850 /im (~ 213.5 /im rest-frame) flux density 
and the assumption of an opacity-weighted blackbody 
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spectral energy distribution (SED) 
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,(i/;Td,/3,z/o) = B,{iy;Ta) (l - exp {-(i^/i^of}) (1) 

Here Td is the dust temperature and i>o is the rest fre- 
quency at which the dust opacity is unity. We will (for 
now) neglect any contribution from hotter dust at wave- 
lengths below the blackbody peak, making the frequency 
integral of this a conservative lower limit on Lir. 
Combined with our upper limit on CO(3~2) emission, 
a given choice of (T^,/?, I'o) then gives us a conservative 
lower limit on the ratio Lin/L'^Q^g ^j- For a particular 
Td, however, this ratio also has a firm upper limit. In 
the molecular interstellar medium of normal galaxies, a 
CO emission line will be optically thick at some bright- 
ness temperature Ti, (less than the gas kinetic tempera- 
ture Tkin as a result of subthermal excitation), while dust 
with Td > Tf, produces far-IR continuum emission that is 
optically thin. As density and column density increase, 
L'qq will rise as the CO transition is thermalized (i.e., 
Tfc Tkin), but Lijt will rise much more dramatically 
as the dust emission becomes optically thick. The upper 
limit on Li^/L'^q is reached when both dust and gas are 
radiat ing as blackbodies with Td = T^. Solomon et al. 
(jl997|) have shown that this "blackbody [upper] limit" is 
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and that local ULIRGs fall only a factor of 3 below it. 

For WMDll, the observed 850 /im flux density and 
upper limit(s) on CO(3-2) emission are only compatible 
with the upper limit on Lir/L^q imposed by Equation 
12 for a severely Hmited range of SEDs. For Avco = 
121 (276)kms-\ dust temperatures T<j < 37 (60) K are 
excluded at once because the blackbody limit they im- 
pose is simply too stringent. More generally, we require 
a quite small value of vq for any choice of Td and f3 if we 
are to suppress Lir, relative to ipo- Figure H shows the 
dependence of c/vq as a function of Td, /?, and Avqo 
if WMDll exactly satisfies the blackbody limit. For 
Avco = 121kms-\ < 60K and (3 > 1.5 will pro- 
duce dust emission that is optically thick at rest wave- 
lengths < 213.5 /xm. Given that local ULIRGs do not 
themselves reach equality in Equation |21 this conclu- 
sion should hold robustly for most plausible {Td, (3). For 
Avco = 276kms~^, dust emission at rest wavelength 
213.5 /im will be optically thick for any SED. 

If WMDll 's observed 850 /xm emission is indeed opti- 
cally thick, we can estimate the characteristic size of the 
emitting region (assuming a spherical geometry) as 
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(70K,1.5,lTHz), WMDll's 
12 T _ and 



(6.0 ± 1.4) X 10^2 /i-^Lq 
380 /i(^7pc. 



Assuming Td = 



For a fiducial (Td, /3, vq 
5*850 translates to Ljr : 
thus to a radius Rd - 

40 and 90 K yields Rd — 610 and 320/i,^ypc respec 
tively (although recall that Td < 60 K is incompatible 
with Avco — 276kms~i). Such radii are compara- 
ble to the sizes of molecula r gas disks in loc al ULIRGs 
(jDownes fc Solomon 1998: .Brvant fc Scoyilldh999}. In- 
deed, the same logic works in reverse: iSolomon et al.l 
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Fig. 1. — Rest wavelength at which WMDll's dust emis- 
sion reaches optical depth unity, assuming the blackbody limit on 
^ir/^CO (Equation 13 is exactly satisfied, as a function of and 
13. Solid (dashed) curves assume Avqq = 121kms-l (276kms-l). 
If WMDll has Lir/Lqq below the blackbody limit, curves will 
shift upwards in this plot. Dot-dashed line shows rest wavelength 
corresponding to observed 850 ^m for z = 2.98. 



(jl997l) point to the compact molecular gas structures in 
ULIRGs as evidence for high optical depths, which they 
use in turn to explain these systems' high Lm/L'^Q ra- 
tios. 

Our conclusion that the dust source in WMDll is both 
optically thick out to long rest wavelengths and spatially 
compact is independent of whether its Ljr is powered by 
star formation or accretion. Given the lack of evidence 
for an active galactic nucleus (AGN) in either optical 
spectroscopy or radio continuum mapping of WMDll 
(jChapman et al.. 2002) . a compact starburst seems more 
likely. However, the cautionary exam ple of the z = 
3.09 Lya nebula detected at 850 //m bv iChapman et all 
(2001), for which accretion remai ns one of the few plau- 
sible sources of ionizing photons ijChapman et al.l l2004'l , 
suggests a deeply embedded AGN should not be excluded 
for WMDll either. 

4. LEGS HAVE DIVERSE MODES OF STAR FORMATION 

Given WMDll's morphological resemblance to an 
early-stage galaxy merger and its large IR lumi- 
nosity, we estimate its molecular gas mass using 
the Mgas/i'co conversion factor ap propriate for lo- 
cal ULIRGs. Ipowncs fc SolomonI (p998) suggest 
that this should be lower than the Galactic value 
and sa 0.8Mq (Kkms"^ pc^)"!, which for Awco - 
121 (276)kms~^ gives WMDll a molecular gas mass of 
< 3.9(5.9) X 10^ Hq^Mq (including helium) that could 
plausibly have orig inated in a single progenitor spiral 
i)Helfer et all \2003f\ . WMDll's ERO component may 
be the remnant of a seco nd, older spheroid progeni- 
tor: IChapman et all ll200S 
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note that if all of its redness 
Ks = 6.15) is attributed to extinction, the sys- 
tem's total 850 /im flux density would be overpredicted 
by a substantial factor. 

To estimate WMDll's star formation rate, we calcu- 
late its total (8 — 1000 /im) Ljr by smoothly grafting 
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a sho rt-wavelength power law Sy oc v ^'^ l)Blain et al.l 
I20n2|) onto the fiducial SED used in P This step prop- 
erly includes emission h'om hotter dust. For Ai>co — 
121kms~^, we begin with = 50K and obta in Ltr = 
(4.0± 0.9) X 10^2 Lq, which for a 1- 100 M^ lSaloeteil 
l^l9Ri^^ Initial Mass Function (IMF) corresponds to a 
star fo rmation rate of (330±80) h^"^ Mq yr~^ ()Kennicuttl 
^23). For Avco — 276kms~^, we must adopt a higher 
Td = 70 K, which gives im = (8.5 ± 2.0) x lO^^ ft-^ Lq 
and a star forma tion rate (710 ± 170) /ig 7 Af© yr~^. 
iPettini et all lj200lD estimate (for the same IMF) that the 
star formation rate in WMDll is 8 h^'^ Mq yr~^ from II/3 
line and rest-UV continuum fluxes before correction of ei- 
ther for extinction. In hindsight, the agreement between 
those estimates does not point to a low global dust ob- 
scuration, but to a complex geometry in which blue star- 
forming patches are scattered about the periphery of a 
dusty starburst core. Such configurations are also seen 
locally: local ULIRGs can have quite blue UV colo rs that 
conce al their true degree of obscuration (Goldade r et al.l 
120021) . while the most massive molecular gas concentra- 
tions and intense star formation can occur off the opti- 
cally bright nu clei in nearby merg ers like the Antennae 
HVigroux et aUllflflfit IWilson et a,].ll200m . 

The ratio of the molecular gas mass and star forma- 
tion rate in WMDll nominally limits the gas exhaustion 
timescale to a short tgas ^10 Myr, modulo large uncer- 
tanties in SED and Avqo- An instructive contrast can 
b e drawn with cB 58, which is forming stars according to 
a lSchmidH l)1959H law with a much longer tgas ~ 240 Myr 



comparable to its most likely age for past star formation 
(B04). Differing modes of star formation in these two 
LBGs naturally explain why we are apparently catch- 
ing WMDll closer to the endpoint of its evolution than 
cB58. In a major merger, gas from the progenitor galax- 
ies will be rapidly dumped into the system's center of 
mass or ejected into tidal tails (from which infall and 
ensuing star formation will be slow). Since the central 
starburst will shut off abruptly as soon as the central gas 
concentration has been consumed or expelled by a wind, 
we might expect that if WMDll is a major merger, we 
should indeed be catching it near the end of a brief burst. 
For more quiescent star formation in a disk, an extended 
H I reservoir would be capable of refilling a molecular 
gas reservoir at smaller radii. We would therefore ex- 
pect that LBGs like cB58, in which star formation has 
a less violent trigger, would tend to be observed in the 
middle of a star formation episode without a sharply de- 
fined end. With this new evidence that LBGs can have 
quite different modes of star formation, it is clear that 
understanding the distribution of such modes across the 
full LBG population will be important in constraining 
these systems' contributions to the cosmic histories of 
star formation and mass assembly. 
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